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Abstract. The purpose of this paper is to prove that the spectrum of the non-self- adjoint one-particle 
Hamiltonian proposed by J. Feinberg and A. Zee (Phys. Rev. E 59 (1999), 6433-6443) has interior 
points. We do this by first recalling that the spectrum of this random operator is the union of the 
set of 1°^ eigenvalues of all infinite matrices with the same structure. We then construct an infinite 
matrix of this structure for which every point of the open unit disk is an eigenvalue, this following 
from the fact that the components of the eigenvector are polynomials in the spectral parameter whose 
non-zero coefficients are zbl's, forming the pattern of an infinite discrete Sierpinski triangle. 
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1 Introduction and Notations 



In this paper we study infinite matrices of the form 



1 
6-1 1 
6o 



\ 



(1) 



with hk G {±1} := { — 1, +1} for ah k G Z. We think of (1) as a hnear operator acting via matrix- 
vector multiphcation on ^^(Z), the standard space of bi-infinite complex sequences with p G [1, oo]. 
By {±1}^ we denote the set of all sequences b = {bk)kez with bk G {±1} for all /c G Z, and we 
refer to the operator on -^^(Z) that is induced by the matrix (1) as A^. For convenience, we will 
also refer to the matrix (1) as A^. For p G [1, oo] and b G {±1}^, we write 



spec^A^ 



spec: 



point 



A' 



{AG 
{AG 
{AG 



A^ - XI is not invertible on ^^(Z)}, 
A^ - XI is not Fredholm on F(Z)}, 
A^ - XI is not injective on F(Z)}. 



Because (1) is a band matrix, it holds (see [23] and [28]) that spec^A^ and specfggA^ do not depend 
on p G [1, oo], so that it makes sense to abbreviate these as spec and spec^gg A^ in what follows. 
Note however that the set of eigenvalues, speCp^-^^ A^, does depend on p. 



Physicists have studied the operator as the (non-self-adjoint) Hamiltonian of a particle 
hopping (asymmetrically) on a 1-dimensional lattice [15, 16, 9, 22]. Applications of such and 
related Hamiltonians, especially examples with random diagonals, include vortex line pinning in 
superconductors and growth models in population biology. The particular model (1) was proposed 
by Feinberg and Zee in [15], and some properties of its spectrum have been studied in [9, 22] (also 
see Paragraph 37, in particular Figure 37.7c, in [38]). 

In all these studies the focus is on the case of a random sequence 6 G {±1}^. A related but 
completely deterministic concept is that of a pseudo-ergodic sequence. In accordance with Davies 
[11], we call b G {±1}^ pseudo-ergodic if every finite pattern of ±l's can be found somewhere (as 
a string of consecutive entries) in b. If b is pseudo-ergodic (which is almost surely the case if all 
G Z, are independent (or at least not fully correlated) samples from a random variable with 
values in {±1} and nonzero probability for both +1 and —1) then, as a consequence of [7] (also 
see [6, 8, 29, 30] and cf. [11]), it holds that 



spec A = speCgggA 



y spec A"" 

ce{±i}^ 



ce{±i}^ 



(2) 



The contribution of [7] is the third sign in (2) that enables, or at least simplifies, the explicit 
computation of the spectra of particular pseudo-ergodic operators in [6, 8, 29]. The first 
sign in (2) follows immediately from the second; the second comes from the Fredholm theory of 
much more general operators and is typically expressed in the language of so-called limit operators 
[34, 35, 27, 8]. (A similar equality, often with the closure taken over the union of spectra, can be 
found in the literature on spectral properties of Schrodinger and more general Jacobi operators 
[32, 4, 10, 11, 21, 1, 31, 17, 18, 19, 20, 33, 26, 25, 36, 37]. The three last papers also shed some 
light on the role of limit operators in the study of the absolutely continuous spectrum.) 

Note that, by (2), the spectrum of A^ does not depend on the actual sequence 6 - as long 
as it is pseudo-ergodic. In [6] we obtain information about the spectrum, pseudospectrum and 
numerical range of the bi-infinite matrix operator A^, its contraction A^ to the positive half axis 
(a semi-infinite matrix) and the finite sections A^ which, for n G N, are n x n submatrices of (1). 
Explicitly and precisely, these related matrices are 



1 



&2 



V 



/ 



/ 



and 





h 



\ 



1 



V 



where in the case n = 1 we set A\ = (0). We explore in some detail in [6] the interrelations between 
the spectra and pseudospectra of A^, and A^. Here, for £ > and a bounded operator A on 
^^(N) or ^^(Z), or on C"^ equipped with the 2-norm, we define the e -pseudospectrum of A (see e.g. 
[2, 38]) by 

spec^A := {A G C : A G spec A or ||(A - A/)"^|| > l/e}, 

where || • || is the induced operator norm. It is convenient also to use the notation specg^ := spec A. 
Note that the finite matrix A^ only depends on the n — 1 values 6i, bn-i G {±1}. Recognising 
this, we will use the notation , where b' = (6i, ...,6^_i) G {±1}^~^, as an alternative notation 
for the same matrix A^. 

Here is a summary of our results from [6]: 



Theorem 1.1 [6] If b G {±1}^ is pseudo-ergodic (which holds almost surely ifb is random in the 
sense discussed above) then the following statements hold. 

a) spec A^ is invariant under reflection about either axis as well as under a 90^ rotation around 
the origin. 

b) Provided the ''positive" part of the sequence b (by which we mean {bk)kef^) itself pseudo- 
ergodic (contains every finite pattern o/ ±1 ^s), then, for all e > one has 

spec^A^ = spec^A^. 
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c) The numerical range of (considered as an operator on £'^{7j)) is 

W{A') = {x^iy : x.yeR, \x\ + \y\ < 2}, 
and spec is a strict subset of the closure, cios {W{A^)), of the numerical range, so that 
spec A^ g {x ^iy : x,y eR, \x\ + \y\ < 2}. 

d) For every n G where {c G {±1}^ : c is n -periodic} , the set 

TTn := U ^P^^^' = U SP^Coint^' (3) 

ceUn ceUn 

is contained in specA^^ by (2). Each set iTn consists of k analytic arcs (see Figure 1.1) 
with ^ < k < 2^ that can be computed explicitly ( as unions of sets of eigenvalues of n x n 
matrices). In particular, 

TTi = [-2, 2] U [-2i, 2i] and 7r2 = tti U {x ^ iy : -1 < x < 1, y = ±x}. 

e) For all n eN and e >0, the set 

an,s := U spec^A^ (4) 

cG{±l} — 1 



IS 



contained in spec^ A^ (see Figure 1.2 for £ = {)). 



f) In the case of spectra (£ = {)), the finite matrix eigenvalues Gn := crn,o from (4) are connected 
with the periodic operator spectra iVn from (3) by 

CTn C 7r2n+2 C SpCC (5) 

for all n (see Figure 2.1). 

g) ^45 a special case of a much more general spectral inclusion result from [5], we can complement 
the inclusion On^e C spec^ from e) by 

o'n,£ C spec^ C (Jn,£+£^ and an C specA^ C clos(<Jn, , 

for n e N and e > 0, where Sn = 4sin6>n < 27r/(n + 1); with On the unique solution in the 

interval ( —r—- — r, , ^ — r ) of the equation 
V2(n + 3)'2(n + l); ^ ^ 

2cos((n + l)i9) = cos((n- l)i9). 

Remark 1.2 Note that the right "=" sign in (3) holds because specA^ = spec^-^^ for ah 
periodic sequences c, whereas the right sign in (2) only holds as stated, with the union taken 
over all c G {±1}^; the spectrum and point spectrum of A'^ are different, in general, for specific 

Remark 1.3 The inclusions in g) imply that specA^ C clos (crn,^^) C clos (spec^^ A^). Since 
6n ^ so that clos (spec^^ A^) specA^ in the Hausdorff metric [38] as n ^ oo, it follows that 
clos(<Jn^£^) specA^ as n ^ oo. For small values of n the upper bound clos(<Jn^£^) can be 
evaluated very explicitly. In particular, = 7r/6 so that £1 = 2 and, since A\ = (0), we obtain 
that specA^ C clos (cri^^ J = {A G C : |A| < 2}. The result in c) above, that spec^d^ is a strict 
subset of the closure of the numerical range, comes from the bound in g) applied with n = 34, 
when cTn^e^ is the union of the pseudospectra of 2^^ ^ 8.6 x 10^ matrices of size 34 x 34. □ 
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Fi gure 1.1: Our figure shows the sets tt^, as defined in (3), for n = 1, ...,30. Recall that tti = [—2,2] U [— 2i, 2i] 
and that, for each n, tti C TTn C {x + iy : x, y G M, \x\ + |y| < 2}. Note also that spectra of periodic infinite matrices 
can be expressed analytically (by Fourier transform techniques, see e.g. [3, 12]) and that each set Tin consists of k 
analytic arcs, where < k <2'^. 

It follows from Theorem 1.1 d) and e) that both 

oo 

(Too '= [J CTn and TToo 

n=l 

are subsets of specA^ (with a^o C tToo, by (5)). These subsets consist of countably many points 
and countably many analytic arcs, respectively, and so both have zero (two-dimensional) Lebesgue 
measure. Indeed, it is not clear from any of the results in Theorem 1.1 (or other results in the 
literature) whether specA^ has positive Lebesgue measure, in particular whether it has interior 
points. Related to this question, Holz et al. [22, Sections I, V, VI], conjecture that clos(croo) C 
spec has a fractal dimension in the range (1, 2), and so has zero Lebesgue measure. 




oo 
n=l 
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The purpose of the current paper is to shed Ught on these questions by constructing a sequence 
c G {±1}^ for which spec^-^^ contains the open unit disk. As a consequence of formula (2) 
and the closedness of spectra, this shows that specA^ contains the closed unit disk and therefore 
has dimension 2 and a positive Lebesgue measure. This is the main result of the next section. 
Intriguingly we will see that the sequence c constructed, while rather irregular, is such that each 
A in the unit disk is an eigenvalue of with an eigenvector u G £^(Z) whose components are 
polynomials in A with coefficients forming the regular self-similar pattern of a discrete Sierpinski 
triangle (7). 

We will finish the paper with our own conjecture on the geometry of clos (ctoo) and spec A^. 



. i . 



• i • 



. I . 

• i ■ 
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Fi gure 1.2: Our figure shows the sets an := <Jn,o of all n x n matrix eigenvalues, as defined in (4), for n = 1, ...,30. 
Note that in the first pictures (with only a few eigenvalues), we have used heavier pixels for the sake of visibility. 
By (5), each of the sets with n = 1,2, 14 in this figure is contained, respectively, in the set number 2n + 2 of 
Figure 1.1. 
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Fi gure 1.3l This is a zoom into (725 ~ the 25th picture of Figure 1.2. The location of this zoom is near the point 
1 + i, which is the midpoint of the northeast edge of the square clos {W{A^)) = conv{2, -2, 2i, -2i}. The picture 
clearly suggests self-similar features of the set (725- 



2 A sequence c for which specA^ contains the unit disk 

The formula (2) for the spectrum of when b G {±1}^ is pseudo-ergodic motivates the fohowing 
approach to decide whether a given point A G C is in spec or not: look for a sequence c G {±1}^ 
such that A G spec^-^^ A*^, in other words, such that there exists a non-zero u G -^^(Z) with 
A'^u = Aii, i.e. 

Ui-^i = Xui — Ci Ui-i (6) 
for z G Z. If such a sequence c exists then A G spec - if not, then not. 
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Fi gure 2.1: Here we see the inclusion (75 C 7ri2, which holds by (5) with n = 5. (The points in (75 are indicated 
by circled dots.) 

Starting from uq = and ui = 1, we will successively use (6) to compute Ui for i = 2,3,... 
(an analogous procedure is possible for i — —1, —2, —3, ...) and see whether the sequence remains 
bounded. Doing so, we get 

U2 = A, ii3 = A^-C2, = A^-(C2+C3)A, 

= A^ - (C2 + C3 + C4)A^ + C2C4, 

and so on. Explicitly, it is easy to check that, for z > 3, the solution of (6) with initial conditions 
ixo = and iii = 1 is given by the characteristic polynomial 

A -1 



-Q-i A 

Thus, for i > 3, is a polynomial of degree i — 1 in A with coefficients depending on C2, 
We will aim to achieve that be a bounded sequence at least for |A| < 1. 
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With this in mind we should try to keep the coefficients of these polynomials small. Precisely, 
our strategy will be to try to choose ci,C2, ... G {±1} such that each Ui is a polynomial in A with 
coefficients in { — 1,0, 1}. The following table, where we abbreviate —1 by ", +1 by and 
by a space, suggests that this seems to be possible. 



i 


Ci 


1 2 3 4 5 ( 


coefficients of A-^ 
) 7 8 9 10 11 


^ in the polynomial Ui 

12 13 14 15 16 ••• 


1 


+ 


+ 






2 


+ 


+ 






3 




+ 






4 




+ 






5 


+ 


- + + 






6 










7 


+ 


+ 


+ 




8 






+ 




9 


+ 


+ 


+ + 




10 
11 


+ 




+ 




12 


+ 


+ - - 


+ 


+ 


13 




+ 


+ 


+ 


14 






+ 


+ 


15 


+ 






+ + 


16 








+ 



(7) 



For i, j G N, denote the coefficient of A-^~^ in the polynomial Ui by pij. Table (7) shows the 
values Pi J for i, j = 1, 16, given the specific choices indicated on the left hand side of the table 
for the coefficients q. From (6) it follows that 

Pi+i,j = Pij-i - CiPi-ij, (8) 

for z G N and j = 1, 2, z + 1, where we have defined pij := if j < 1, z < 1, or j > i. 

Let us explore more systematically whether it is possible to choose the coefficients q so as to 
ensure that all the coefficients pij G { — 1, 0, 1}. Note first that, if this is possible, then if, for some 
i, j, one has that Pij-i 7^ and Pi-ij 7^ 0, then Pij-i^ Pi-iJ ^ {~^^ Thus it follows from (8) 
that Pi-\-ij = 0, i.e. 

Q = Pij-i/Pi-ij = Pij-iPi-ij, (9) 

since otherwise Pi-\-ij G {—2,2}. Illustrating this, look at pi5^i = —1 and pi4,2 = —1 in the above 
table. If we chose C15 = —1, we would get from (8) that pi6,2 = —2 { — 1,0, 1}, so it is necessary 
to choose Ci5 = 1 = pi5,i Pi4,2- Luckily, the same value C15 = 1 is required by the values of ^15^9 and 
Pi4,io, as well as by ^15,13 and ^14,14. We will prove that this coincidence, i.e. that the right-hand 
side of (9) is (if non-zero) independent of j, is not a matter of fortune. As a result we will show 
that the pattern of coefficients in table (7) continues without end, only using values from { — 1, 0, 1} 
for pij and from {±1} for q. To prove this, we will make use of a particular self-similarity in the 
triangular pattern of (7); more precisely, we will show that the pattern of non-zero values of the 
coefficients pij forms a so-called infinite discrete Sierpinski triangle. 

Proposition 2.1 Define the sequence c G {±1}^^ for positive indices by ci = 1 and by the require- 
ment that 

C2i = C2i-iQ and C2i+i = -C2i, z = l,2, 
and for non-positive indices by 

c-i = Q+i, i = 0, 1, ... . 
Further, given A G define the sequence u = {ui)i^z, by the requirement that 

Ui^i = Xui - CiUi-i, i G Z, 

and by the initial conditions 

Uq =0, Ui = 1. 
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Then, as a function of X, for i G Ui is a polynomial of degree \i\ — 1 with all its coefficients 
taking values in the set { — 1,0, 1}. 

In more detail, denoting, for i^j G the coefficient of \^~^ in the polynomial Ui by Pij, the 
following statements hold. 

('^) Pi J = /^^ j > h ^0 that, for every i G 



Ui = ^PijX^ ^ 

(a) Defining, additionally, pij := z, j G N U {0} and i = or j = 0, it holds that pi^i = 1 
and that 

Pi+ij = Pij-i - CiPi-ij, (10) 

for i and j = 1, 2, i + 1. 
(Hi) pij =0 if i -\- j is odd. 

(iv) Writing the semi-infinite coefficient matrix P = {pij)ijeN i'^ block form as P = {pij)ij^fq 
where, for i^j G N; 

T). . '= ( ^2i-l,2j-l P2i-l,2j 
''^ ' V ^2i,2j-l P2i,2j 

it holds, for i gN, that Pij = for j > i while, for j = 1, 



(v) Pij G {-1,0,1} forij G N. 



1 

P^.[ 1 

1 

C2i-lPi-lJ [ Q Q 



if i -\- j is even^ 
if i -\- j is odd. 



(11) 



(vi) Let V := {(0,0), (-1,-1), (1,-1)} and let S := {{ij) G : Pij is non-zero}. Let S := 2^' 
be the set of all subsets ofN'^, and define T : H ^ E by 



Then, where 
it holds that 



T(cr) := 2a + V = {2a ^ b : a e a, b e V}, forae^. 

Si := {(1, 1)}, and S^+i := T(^^), n G N, 
5' = y 5'n and that S = T{S). 



nef' 



(vii) Fori G N U {0}, 

where, for j G N U {0}, 



di Ui, 



2j 



2i 



@\ 


\ i 




X 






X 






y 





Remark 2.2 Statements {iv) and {vi) reveal the self-similar nature of 
the pattern (7). With respect to a scaling of the pattern by the factor 2, 
an entry pij, with i -\- j even, replicates three times: as P2i-i,2j-i^ P2i,2j 
and, multiplied by C2i+i, as P2i+i,2j-i- So the "volume" of the pattern (7) 
triples under a scaling by 2, which is why (see [13]) its zeta dimension is 
log2 3 ^ 1.585 - exactly the fractal (Hausdorff or box-counting) dimension 
of its bounded version, the usual Sierpinski triangle or gasket [14]. □ 



As an immediate consequence of Proposition 2.1 and formula (2) we get our main result. 
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Theorem 2.3 For the sequence c G {±1}^ from Proposition 2.1, it holds that the closed unit disk 
B := G C : |z| < 1} is contained in spec . Consequently, for a pseudo-ergodic b G {±1}^^ one 
has D C specA^, so that spec has dimension 2 and a positive Lebesgue measure. 



Proof. Let A G B := G C : \z\ < 1}, let c be the sequence from Proposition 2.1 and u : Z ^ 
the corresponding eigenfunction from (6). Then, for every z G N, 



'I oo ^ 



since pij G { — 1, 0, 1} for ah z, j, showing that u G ^^(Z), and, by our construction (6), A^u = Xu. 
So D C spec^-j^^ C spec A^. Since spec is closed, it holds that D C spec A^. The claim for 
a pseudo-ergodic b now follows from specA^ C specA^, by (2). Finally, from the monotonicity of 
(all notions of) dimension [14], it follows that 2 = dim(D) < dim(spec A^) < dim(R^) = 2. ■ 

Proof of Proposition 2.1. Statements {i) and (n) are clear from the discussion preceding 
Proposition 2.1, and statement (Hi) then follows easily by induction. Thus Pij = for j > i, and 
in every matrix pij the off-diagonal entries are zero, i.e. P2i-i,2j = = P2i,2j-i for all j. 

We will now prove {iv) by proving by induction that, for each z G N, (11) holds for j = 1, 
It is easy to check that (11) holds for z = j = 1. Now suppose that, for some k G N, (11) holds 
for i = l,...,/c, j = l,..,i. We will show that this implies that (11) holds for i = /c + 1 and 
j = l,...,/c + l. 

We let i = /c + 1 and start with the case when i -\- j is even. By (10) we have that 

P2i-l,2j-l = P2i-2,2j-2 — C2i-2P2i-3,2j-l 

= P2(i-l),2(j-l) - C2i-2P2(i-l)-l,2j-l, (12) 

with P2(i-i),2(j-i) = Pi-i,j-i = if j = 1 and, by the inductive hypothesis (and since i — 1 -\- j — 1 
is even), P2(i-i),2(j-i) = Pi-ij-i if j > 1- Also, by the inductive hypothesis, P2(i-i)-i,2j-i = 
C2{i-i)-iPi-2,j since i — l-\-j is odd, while, from the definition of the sequence c, C2i-2 = C2i-3 
Inserting these results in (12), we get that 

P2i-l,2j-l = Pi-l,j-l — C2i-3 Ci-i C2i-3Pi-2,j 

= Pi-lJ-1 - Ci-iPi-2,j =PiJ^ (13) 

by (10). We have observed already that P2i-i,2j = = P2i,2j-i for all z, j, so it remains to consider 
P2i,2j' By (10), (13), and the inductive hypothesis which implies that P2i-2,2j = P2{i-i),2j = as 
i — 1 + j is odd, we have that 

P2i,2j = P2i-l,2j-l — C2i-lP2i-2,2j = PiJ- 

Now suppose i -\- j is odd. Then, by (10) and the inductive hypothesis, 

P2i-l,2j-l = P2i-2,2j-2 — C2i-2P2i-3,2j-l 

= - C2i-3Ci-iPi-ij = C2i-lPi-lJ, 

since C2i-i = — C2i-2 = —C2i-3 Ci-i. By (10) and the inductive hypothesis and noting that i — l-\-j 
is even, 

P2i,2j = P2i-l,2j-l — C2i-lP2i-2,2j 

= C2i-lPi-lJ - C2i-lPi-lJ = 0. 

This completes the proof of {iv)^ and {v) follows from {iv) by a simple induction argument. 

To see that {vi) is true, observe first that, from (i), (in), and (iv) (and cf. Remark 2.2), it 
holds for i'J' G N that {i' J') G S iff, for some iJeN either {i' J ) = (2i,2j) and {ij) G S; or 
(i',f) = {2i - l,2j - 1) and (ij) G or {i\f) = {2i + l,2j - 1) and (ij) G S. From this it 
follows that 5' = T{S). 
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Define a metric d on S by 

(i,j)G(aUr)\(anr) 

Tiien, since (T(cr) U T(r)) \ (T(cr) fl T(r)) C T {{a U r) \ {a H r)) for all a, r G S, 

(i,j)G(aUr)\(anr) 

2-'-^ {2-'-^ + 2^-^-^' + 2-'-^) < ^ d{a, r), (14) 

(i,j)G(crUr)\(crnr) 

if (1, 1) (cr U r) \ (cr n r). Let Si := {a G S : (1, 1) G a}. Then T(Si) C Si and, by (14), T is 
a contraction mapping on Hi. Thus, by the contraction mapping theorem, T has a unique fixed 
point in Si, which is the set 5, and, if <Ji G Si and a^+i := T(<Jn), n G N, then d{(Jn^S) 
as n ^ oo. In particular, d{Sn^S) ^ as n ^ oo. Since also (by an easy induction argument) 

6*1 C ^2 C it follows that S = UnGN^n- 

Define V-i for z = 0, 1, ... hy V-i := diUi, which implies that Vq = 0, and set ^i = 1. Then, since 
Ui is defined uniquely for i < by the requirement that it satisfy (6) for z < with the initial 
conditions that uq = and iii = 1, to show {vii) it is enough to check that the sequence Vi satisfies 
(6) for z < 0, i.e. that 

v-i+i = Xv-i - C-iV-i-i, ^ = 0, 1, ... . 
But vi — Xvq + cqV-i = 1 + codiUi = 0, so the equation holds for z = 0, and, for i G N, 

v-i^i - Xv-i + C-iV-i-i = di-iUi-i - XdiUi + Q+i^i+iiii+i 

= {di-i - CiCiJ^idiJ^i)Ui-i - X{di - Ci^idi^i)ui, 

since iz^+i = Xui — CiUi-i. Since uq = 0, the right hand side of this last equation is zero for i G N 
provided that di = q+iG^^+i for z G N. But this follows from the definitions of the sequences c and 
d. m 

Remark 2.4 The standard infinite discrete Sierpinski triangle (e.g. [24]) is the set 5 C defined 
by S := UnGN-Sn, where Si := {(1,1)} and the sets Sn^ n = 2,3,..., are defined recursively by 
Sn-\-i '= 2Sn + V, where V := {(0, 0), (— 1, — 1), (0, — 1)}. One instance where S arises is as the 
pattern of odd coefficients in Pascal's triangle: for i G N and j = l,...,z the coefficient of x^~^ 
in (1 + xy~^ is odd iff (i,j) G so that the discrete Sierpinski triangle is often referred to as 
Pascal's triangle modulo 2 (e.g. [13]). Proposition 2.1{vi) (cf. Remark 2.2) makes clear that the 
pattern 5 C of the non-zero coefficients in table (7) is essentially that of the standard discrete 
Sierpinski triangle S; indeed, the sets S and S are connected by a linear mapping: (i,j) G iff 
{2i- jj)eS, for ij eN. □ 

Remark 2.5 Note that the sequence c from Proposition 2.1 is not pseudo-ergodic since, by 

C2i+i = — C2i, the patterns "+ + +" and " " can never occur as consecutive entries in the 

sequence c. □ 

Based on Theorems 1.1 and 2.3 and the numerical results displayed in Figures 1.1 and 1.2, we 
make the following conjecture. 

Conjecture. We conjecture that clos(<Joo) = clos(7roo) = specA^, and that spec is a simply 
connected set which is the closure of its interior and which has a fractal boundary. 
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